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Abstract. In this paper, we suggest a new way to identify single bright sources of 
Ultra High Energy Cosmic Rays (UHECR) on the sky, on top of background. We look 
for doublets of events at the highest energies, E > 6 ■ 10 19 eV, and identify low energy 
tails, which are deflected by the Galactic Magnetic Field (GMF). For the sources which 
are detected, we can recover their angular positions on the sky within one degree from 
the real ones in 68% of cases. The reconstruction of the deflection power of the regular 
GMF is strongly affected by the value of the turbulent GMF. For typical values of 4/iG 
near the Earth, one can reconstruct the deflection power with 25% precision in 68% of 
cases. 
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1. Introduction 

The search for the Ultra-High Energy Cosmic Ray (UHECR) sources is one of the most 
interesting problems of modern astroparticle physics. The observation of the Greisen- 
Zatsepin-Kuzmin (GZK) cutoff p] in the energy spectrum by HiRes experiment [2] 
proves the astrophysical origin of UHECR. 

Several classes of astrophysical sources can potentially accelerate cosmic rays to 
the highest energies [3], but some of them are restricted by different energy losses in 
the sources (for a recent review see, for example, Ref. [3]). Also, for each class of 
sources, several alternative acceleration mechanisms were suggested, including standard 
first order Fermi acceleration in shocks or acceleration in induced electric fields in the 
polar regions of supermassive black holes [5]. 

Accelerated cosmic rays lose their energy during the propagation from the sources 
to the Earth due to pion production, pair production and redshifting. These energy 
losses define the horizon - the typical distance from which cosmic rays can reach us with 
a final energy bigger than some value. The horizon value was discussed in Ref. [6], in 
the constant energy loss approximation. It was updated to realistic stochastic energy 
losses in Ref. [?]. Finally, the effects of the experimental energy resolution were taken 
into account in [8]. Both for protons and iron nuclei, 70% of sources would be located 
within 70 Mpc of the Earth for energies E > 10 20 eV, and within 250 Mpc of the Earth 
for E > 6 • 10 19 eV. 

Since the Large Scale Structure (LSS) is still anisotropic on 100-200 Mpc scales, any 
class of astrophysical sources of UHECR should give an anisotropic signal in detectors 
with typical angular scales of tens of degrees. Anisotropies on medium scales, 20-30 
degrees, have been found previously by combining all the available data of "old" cosmic 
ray experiments [9j. Such anisotropies were predicted as a consequence of the observed 
LSS of matter [10] and favor therefore, with the presence of the Greisen-Zatsepin-Kuzmin 
cutoff, the extragalactic origin of UHECRs. 

These results still did not allow to find the sources of UHECR, both due to limited 
statistics and to the fact that any different class of astrophysical sources is located in 
the same LSS, which makes their identification more difficult. 

Another important factor which complicates the detection of cosmic ray sources is 
the deflection of UHECR in the galactic and extragalactic magnetic fields. Extragalactic 
fields are unknown outside of galaxy clusters. The numerical simulations done by several 
groups predict a wide range of deflections for 10 20 eV protons. According to Ref. [TTj, 
deflections are negligible (5 < 1° in 98% of the sky within 100 Mpc of the Earth). 
On the contrary, according to Ref. [12], deflections are more than 10 degrees in 70% 
of the sky. In this case, the search for the UHECR sources is extremely difficult and 
close to impossible, even for proton primaries. In the following, we will assume that 
the deflections in the extragalactic magnetic fields are small, following the results of 
Ref. [XT]. Let us note that if the sources are located in regions with high values of 
magnetic fields, like galaxy clusters, both their spectrum and their detected flux can be 



Search for single sources of ultra high energy cosmic rays on the sky 



3 



affected by magnetic lensing effects in the cluster [13]. Below, we assume that one can 
neglect such effects. 

Even in case the extragalactic magnetic fields are negligible, UHECR would still be 
deflected in the Galactic Magnetic Field (GMF). For the GMF, the situation is much 
better due to existing measurements of both the regular and the turbulent components 
of this field. However, there is no theoretical model of the GMF which can fit all existing 
data. The measured values of the GMF still allow to estimate the typical deflections 
of UHECR in it, which are ~ 0.5-2 degrees for protons of energies equal to 10 20 eV, 
depending on the position on the sky and on the GMF model. 

There are two possible ways to continue the search for the UHECR sources in the 
near future. One is to look at the highest energies only, E > 6 ■ 10 19 eV, and collect 
enough statistics to find point sources. This method can be strongly affected by the 
presence of a large fraction of nuclei in the UHECR flux. 

In this paper, we present another way to search for individual bright sources. One 
can hope that the regular GMF gives the major contribution to UHECR deflections. 
Then, one can try to look near high energy events for low energy tails on top of a high 
background of UHECR, which comes from all unresolved proton sources and from heavy 
nuclei. In the present paper, we suggest a new method which is based on searching for 
low energy tails near doublets or clusters of events with high energies, E > 10 19,8 eV 
~ 6.3 ■ 10 19 eV. This allows us to reduce the background significantly, compared to 
usual searches for tails near every single high energy event. 

We simulated the signals from a single source and from a background of UHECR. 
Following the results of HiRes [H], we take a proton source, except in a paragraph 
of Section 0], where we investigate some possible extensions of this work to sources of 
heavier nuclei. 

Then, we studied how the probability to detect this source depends on unknown 
parameters of the source and of the magnetic field. In case the source is detected, we 
finally reconstruct its position on the sky and the local deflection power of the regular 
GMF. 

In Section [21 we simulate the signal from a single source on top of background. 
Section [3] depicts the method used to detect the source and reconstruct its position. 
In Section HJ we test how its efficiency varies with experimental or unknown physical 
parameters. In Section [5j we summarize our results. 

2. Simulation of the signal from a single source on top of background 

For this study, we only take into account the cosmic rays which energies are bigger 
than a given energy threshold, E^. Typically, E^ > 10 19 eV. Detecting a signal from 
a source on top of background at lower energies is hopeless for realistic magnetic field 
parameters. 

We start below with the simulation of the source signal. 
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Figure 1. Flux (left) and arrival directions (right) of UHECR protons after 
propagation on 50 Mpc and deflection in the Galactic Magnetic Field. Left: solid 
red line for the power law index a = 2.2 and green dashed line for a — 2.7. Both 
spectra normalized to the same flux above 10 198 eV. Right: projections of the arrival 
directions of cosmic rays emitted by a source (red cross) in the plane tangent to the 
celestial sphere and centered on the highest energy event (blue filled square in (0,0)). 
The deflections X and Y are measured on two orthogonal axes, in degrees. Open 
symbols for the CR from the source, deflected by the regular GMF only. Filled symbols 
for the same CR when the turbulent component is also taken into account. The shapes 
correspond to the CR energies (triangles: below 10 193 eV, circles: between 10 19 3 and 
10 19 - 8 eV, squares: above 10 198 eV). 



2.1. Events from the source 

We assume that the source has a power law acceleration spectrum in the form: 

f(e) = /-Le (E ma , - E) , (1) 

where / is the flux normalization factor, a the power law index, O the Heaviside step 
function and E max = 10 21 eV the maximum energy to which the source can accelerate 
cosmic rays. Below, we consider two examples of power law spectra: a = 2.2 for Fermi 
acceleration and a = 2.7 for V. Berezinsky et al. model [TS"] . 

Since we would like to study a bright nearby proton source, we take typical distances 
of 50-100 Mpc and propagate the protons to the Earth. UHECR protons lose energy on 
their way to the Earth, due to pion production (GZK effect), e ± pair production and the 
expansion of the Universe. We use the results of Ref. [8] for the energy losses. Typical 
spectra are presented in Fig. [Th. for a source located at 50 Mpc from the Earth. The two 
curves are normalized to the same luminosity above E = 10 19,8 eV ~ 60 EeV. For such 
small distances, the propagation affects the injection spectrum only at high energies, 
E > 10 19 ' 6 eV ~ 40 EeV, creating the "bump" [16] and the cutoff in the spectrum. 

The number of events from the source is computed according to its luminosity. It 
is defined as the fraction of source events with energies above 10 19 ' 8 eV in all the sky. 
In practice, the number of source events which are generated above E t h in every Monte 
Carlo simulation fluctuates around the mean value inferred from the source luminosity. 
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On the way from the source to the Earth, cosmic rays are deflected in the Galactic 
Magnetic Field. This field consists of a regular and a turbulent component which are 
different in the disk and in the halo. At high energies, E > 40 EeV, and outside of the 
galactic plane, the deflections of cosmic rays in the regular field S reg , can be estimated 
as (see e.g. Ref. [T7]): 



u reg — 



,40EeV 
E/Z 



L ds B 



o 3kpc 2fiG 



(2) 



In other words, 



$reg — g j (3) 

where T> is the deflection power of the regular field for a given region of the sky. 

The regular component shifts the cosmic rays as 1/E along a curve. This curve can 
be approximately straight outside of the galactic plane for high enough energies. Since 
the magnetic field in the disk is poorly known, one has to exclude the galactic plane for 
single source searches in any case. 

As an example of regular field, we take a model based on the model of Prouza and 
Smida (PS) [i~8|[T9]. The PS model consists of a bispiral field in the disk and a halo field 
containing a poloidal and a torroidal component. We updated this model, according to 
the latest results of Refs. [20TI24"] . and a discussion during the Ringberg workshop [25J. 
For instance, we choose an exponentially decaying profile along the galactocentric radius 
r, instead of a 1/r profile and we set the value of the regular field close to the Sun, B reg , 
to 2 fiG. 

For the turbulent component, we take the model detailed in Refs. [26, 27J. For 
particles travelling in a turbulent field, the distribution of angular deflections is a two- 
dimensional Gaussian centered on zero (no deflection). 

The root mean square deflection due to the turbulent field, S rms , can be written as: 



a/2 E ^ ^ E/Z J \^GJ V 3kpcV 50pc V; 

where B rms is the RMS value of the turbulent magnetic field, L c its correlation length, 
and L ^> L c the length travelled by the cosmic rays in it. We take L c = 50 pc, which is 
a commonly admitted value [25] . 

In this work, we take the following profile for B rms : 

B rms = B(r) exp > ( 5 ) 

where r is the galactocentric radius, and z the distance to the galactic plane. We take 
here Zq — 1.5 kpc. 

The radial profile of the field can be defined as: 

B rms0 ■ exp (|f) , if r < 3kpc (bulge) 

B(r) = { ^ /:(r-8.5kpc)\ . f 01 (6) 



„ /-(r-8.5kpc)\ 01 

B rms0 ■ exp I v g5kpc J , if r > 3kpc 
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where B rms0 denotes the value of B rms close to the Sun. The currently admitted value 
is B rms0 ~ 4 fiG [25]. We take it as a reference value, but we will also discuss about the 
impact of weaker and stronger fields. 

One can note that in a given region of the sky, 5 reg and 5 rms are approximately 
proportional to 1/E. Therefore, for a given energy one has approximately 5 rms oc S reg : 

firms — T^tur-b/reg ' 0~reg (7) 

IZturb/reg estimates the ratio of deflections respectively due to the turbulent and 
regular fields. 

In terms of deflections, the impact of the two GMF components can be seen in 
Fig. [lb. The regular component shifts the cosmic rays as 1/E along a curve. For each 
energy, the turbulent component spreads the cosmic rays around the position where they 
would have been deflected in the regular field only. Due to Eqn. (J7|), the cosmic rays are 
roughly spread in a sector. When there are many events with the same energies, they 
are not spread uniformly, as discussed in Refs. [TTII2S1I2HII29J . Note that our study does 
not depend on the details of the Galactic Magnetic Field model. The results mostly 
depend on two parameters that summarize the impact of the regular and turbulent fields. 
The first of them is the deflection power of the regular field in the region surrounding 
the source, V -see Eqn. ([3l). The second one is the proportionality factor lZ tur b/reg in 
Eqn. ©. 

One can think that the source detection becomes easier when lowering the energy 
threshold E^, thanks to the increasing number of events. However, it is not so, because 
the background grows much faster at low energies. We will discuss this in the next 
section. 



2.2. Background 

We generate the background according to the exposure of Telescope Array, following 
Ref. [HO]- We set the declination of Telescope Array to 5 = 40.1° and its maximum 
zenith angle to 9 = 45°. The energies of background events are generated according to 
the spectrum measured by HiRes [2j. 

Fig. [2] shows source and background events in the region around the source for 
different energy thresholds. In Fig. [2^ we plot the arrival directions of cosmic rays 
with energies bigger than E t h = 10 190 eV, and in Fig. [2b, with energies bigger than 
E t h = 10 19 ' 6 eV. One can clearly notice two important facts for source detection. At 
high energies, on the one hand, we lose in terms of statistics because we have less cosmic 
rays from the source, but on the other hand, the background is more strongly reduced 
than the source signal. This is due to the difference of slopes between the background 
and the source spectra. That is why there should be an optimum energy threshold E t h 
for source detection between 10 19 eV and the very highest energies. 

For deflections lower than 10 degrees, when one can neglect the dependence on 
exposure on the considered portion of the sky, our results do not depend on the position 
of the source because both the signals from the source and the background follow the 
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Figure 2. Projections of the arrival directions of CR in the tangent plane as in 
Fig. [TJd. Red cross for the source position, blue filled circles for source events and 
green open squares for background events. Left panel (a): energy threshold equal to 
Ey, = 10 19 eV; Right panel (b): same picture, but with a higher energy threshold, 
equal to E th = 10 19 6 eV. 



same exposure. The regions of very small exposure are particular cases. We need at 
least 3 events from the source for our analysis and in these regions this condition can be 
difficult to hold for not very bright sources. For the results we discuss in this paper, we 
put the source in a region where the exposure of Telescope Array is not maximal but 
still significant. 

In Section HI we also investigate separately how our results are affected by the 
anisotropy at energies above 10 19 ' 8 eV. 

3. Method to find the source on top of background 

In the previous Section, we described how we generate the background events and those 
coming from the source. Now, we mix all of them in the same "data" file and try to 
find the source on top of the background. For this purpose, we develop a new source 
detection method in the next section. When the source is detected, we reconstruct 
its position on the sky and the deflection power of the regular magnetic field near the 
source, V. 

3.1. Source detection and background rejection 

Usually, one starts with a single high energy event with an energy Ei satisfying, for 
example: Ei > 10 19 ' 8 eV ~ 6.3 • 10 19 eV. This guarantees two things. First, the source 
is located nearby in the local Universe, at distances smaller than 250 Mpc in 70% of 
cases -see Ref. [8]. The source density within such distances should not be too large 
to prevent the identification of the brightest sources. Second, for standard power law 
injection spectra with indexes a > 2, one should see some low energy events from the 
source near this high energy event. 



Search for single sources of ultra high energy cosmic rays on the sky 



8 




Figure 3. Example of source detection on top of background. The big black circle 
with a radius according to the Eqn. ([5]) is centered on the highest energy event, in 
(0,0). The black arrow points at the second highest energy event. The sector in the 
direction of this event is highlighted in grey. Blue circles indicate source events. Green 
squares show positions of background events. Events that pass the final selection 
are surrounded with magenta circles. (X,Y) and (X',Y') respectively denote old and 
new systems of coordinates. S and S' for real and reconstructed source positions, 
respectively. 



Let us consider the region around one high energy event with E x > 10 19 ' 8 eV. Its 
coordinates are redefined as (0,0) in Fig. [3j First of all, we restrict our analysis to 
the cosmic rays which angular distance to this event is lower than a given value. This 
maximal angular distance is denoted by the big black circle in Fig. [3j Its radius R 
depends on the energy threshold E^ and has to be optimized for each E^ and Ei. In 
case the cosmic rays produced by the source are roughly deflected along a line with 
deflections proportional to 1 /E, one can expect that R varies in the following way: 

R(M = R(l0 19 - 8 eV)^-^' ( 8 ) 

where R (10 19,8 eV) = 3—6° depending on the position on the sky, the model of GMF 
and the composition of primaries. This linear dependence is satisfied for approximately 
E t /j > 10 19 6 eV. For lower energies, optimum radii are smaller than radii given by this 
equation. The optimization of all parameters is discussed in Section 13.21 
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If the turbulent field strength B rTOs0 is not too big, most of the lower energy source 
events are deflected in a region that is sector-shaped (see, for example, Fig. [2^). A good 
way to find the direction of this sector is to look for a second highest energy CR around 
the highest energy one. The line containing both events should generally point at the 
tail of lower energy events from the source. Let us denote E 2 (E 2 < Ei), the energy of 
the second highest energy event. 

In order to discriminate between cases when the highest energy event comes from 
the source or from the background, one has to set two requirements: 

(i) Its energy E 2 has to be above a given threshold E 2min : E 2 > E 2min . We will show 
below that when the contribution of the turbulent field is not negligible, one should 
set E 2min to 10 19,8 eV. Otherwise, one is too often confused by detecting fake sources 
in the background. 

(ii) Its angular distance to the highest energy event rf 12 should not exceed a maximum 
value. We use the following condition: 

V V 

da < V " iT (9) 

where f3 ~ 1 is a parameter that is optimized (see below) . This second requirement 
guarantees the compatibility with an emission from a same source. If the angular 
distance is too big, the second highest energy cosmic ray probably comes from the 
background and one cannot trust it to find the direction of the sector. 

Both above conditions are optimized to increase the detection of the source and 
reduce the background of fake sources. If more than one event in the circle fulfils these 
two conditions, then the different potential second highest energy events are tested by 
decreasing energy order. 

The next step is to define the sector. In the following study, we only work with 
the events located in the sector, assuming that the regular magnetic field dominates 
deflections. Its central axis is defined as the line which contains the two highest energy 
cosmic rays. The optimum value of the sector angle mostly depends on the turbulent 
field strength. In Fig. [3j the sector is highlighted in grey, its central axis is the central 
black line and the black arrow points at the second highest energy event. The use of the 
sector enables us to reduce significantly the background. In Section HI we take 9 = 60°, 
which reduces the background by a factor 6. 

Then, we rotate from the old coordinates (X,Y) to new ones (X',Y') -cf. Fig. [3j 
The direction of X' is given by the center of mass of all points in the sector, which is 
not exactly the direction of the central axis of the sector. 

The correlation coefficient between X' and 1/E, Corr(X',l/E), can be defined as 
following [31] : 

0«r(XM/E)= , C0V(X '- 1/E) (10) 
/Var(X')Var(l/E) 
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where Cov and Var respectively denote covariance and variance: 



Cov(X',l/E) = -X;(X' i 



1 ^ 



(X')) (1/E, - (1/EJ) 



(11) 



i=l 



Var (X') = (X' 2 ) - (X') 2 



1 n 

-E(X',-(X)) 2 



(12) 



if n is the total number of points for which the correlation coefficient is computed. 

The next step consists in discriminating the considered selected regions according 
to the value of the correlation coefficient Corr(X',l/E) defined in Eqn. ffTOl . A minimal 
value of this coefficient below which the regions are not detected, C m i n , can be found 
by optimizing the rejection of the background as discussed below. 

The use of the correlation coefficient between X' and 1/E was first suggested in 
Ref. for deflections in the regular field (R-turb/reg "C 1). When the contribution 
of the turbulent field is added, there are two cases: for values of IZturb/reg lower than 
~ 0.25, one can still benefit from the correlation coefficient method depicted below, but 
when TZturb/reg ~ 0.4, Corr(X',l/E) for source events drops to values comparable to 
those for background events (see Fig. [9]). 

The next step is to erase in the considered sector the points which decrease 
Corr(X',l/E). In principle, most of these odd points come from the background. The 
CR that are located in the sector are tested one after another, by increasing energy 
order. We do so because lower energy CR more probably come from the background. 
Each tested CR is temporarily removed. Then the new axes (X',Y') are redetermined, 
because the direction of X' is defined by the center of mass of the considered points in the 
sector, and Corr(X',l/E) is recomputed. If the new value of Corr(X',l/E) is lower than 
the previous one, the tested point is rewritten in the list of CR and is definitely kept. 
Otherwise, it is definitely removed. In Fig. [3j the events that are kept are surrounded 
with magenta circles. These tests stop in two cases. First, if Corr(X',l/E) > C m j n . 
Second, if all the events in the sector have been tested, or if there are only 3 points left 
because Corr(X',l/E) would not make any sense with less than 3 points. 

If the final correlation coefficient is larger or equal to C m j n , the considered selected 
region around the highest energy event is detected: if it was a source event, the source is 
detected. If Corr(X',l/E) < C m j n , and if there is another second highest energy event, 
the program tries the new associated sector and redoes the depicted method. Otherwise, 
the region is rejected and there is no detection. 

We tested if we would have better results if we only check the points at the border 
of the sector, those which are far from the central axis, instead all of them. But if O is 
correctly optimized, no significant difference could be seen on the results. 

3.2. Optimization of the parameters of the method 

In the previous section, we defined five parameters which must be optimized: R, C min , O, 
j3 and ^min- They are "internal" parameters of the method. Their optimum values can 
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Figure 4. Dependence on the minimum energy for the second highest energy event, 
E>2min- Green dotted and black dashed lines for the probabilities to have at least one 
second highest energy event above "&2min around a source or a background highest 
energy event, respectively. Red thick and blue thin solid lines for the probabilities 
to detect respectively the source and the background through the procedure. Source 
located in a region of the sky where T> ~ 5.3° x 40 EeV and Ti-turb/reg — 0.2. Source 
of luminosity equal to 2.6%, Eth — 10 19 3 eV, 5000 events above 10 19 eV in the whole 
visible sky, B re9 = 2 /iG and B rms o = 4/iG. 

depend on physical or experimental parameters, as the regular and turbulent magnetic 
fields strengths, the energy threshold of sky maps, etc. Optimization is not unique. It is 
a compromise between the acceptance of the signal and the reduction of the background. 
Besides, the optimized parameters are not independent. 

The optimization of E 2m j n , the minimum energy for the second highest energy event, 
is shown in Fig. HI When the turbulent field strength is not negligible, one has to require 
that E 2min = 10 19 ' 8 eV. If we set E 2min to lower values, the probability to detect the 
source is overwhelmed by the probability to detect some background. The detection of 
the source tends to become hopeless for E 2m ; n below 10 19 ' 8 eV. Therefore, we will only 
search for at least doublets of events above 10 19 ' 8 eV in the sky. In the following, the 
word "doublets" (from the source or from the background) will implicitly refer to such 
doublets of events with energies above 10 19 ' 8 eV. 

If the turbulent field is negligible as in Ref. [29], one could take a smaller sector 
angle 6 and use a lower value of E 2min . 

In Fig. [5j we give examples of plots used to optimize the four other parameters. 
They are done for a source which luminosity is equal to 2.6% of the total luminosity in 
UHECR above 10 19,8 eV. It corresponds to the absolute luminosity L ~ 1.5-10 41 erg-s -1 
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Figure 5. Dependence on the four other internal parameters, at Ey, = 10 19 6 eV. 
Top left corner (a): radius of the circle; top right (b): coefficient C m i n ; bottom 
left (c): sector angle 0; bottom right (d): coefficient j3. Magenta dashed-dotted 
line and red thick solid line for the probabilities to have a source doublet and to also 
detect it, respectively. Black dashed line and blue thin solid line for the probabilities 
to have a background doublet and to also detect it, respectively. Green dotted line 
for source doublets and clusters. Difference between doublets and clusters explained 
in Section |U 



for E > 10 19 ' 8 eV and for a source located at 100 Mpc from the Earth. We used, as an 
example, the measured flux and the detector parameters of Auger experiment [32]. We 
also fixed here E t h = 10 19 ' 6 eV and B rms o — 4 /xG. For this particular example, there 
is a range of parameters in which the values are close to optimal. 

In Fig. El we can take for example j3 = 1.25, because the probability to have a 
doublet from the source is already maximal for this value. Taking a higher value would 
just increase the background. The best value for the opening angle of the sector is in 
the range 6 ~ 60° - 80°. 

3.3. Reconstruction of the source position 

When the source is detected with the method discussed above, the program reconstructs 
its position, and computes the local deflection power of the regular field in the 
surrounding region, T>. To do so, we plot 1/E versus X' for all the events which are 
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recognized as source events and fit it with a straight line. From this fit, we can compute 
the deflection power and reconstruct on the axis X' the position of the source -see Fig. |3j 

4. Results 

In this Section, the results are presented for a source which is located in a region where 
V ~ 5.3° x 10 19 - 6 eV ~ 5.3° x 40 EeV and n turb/reg ~ 0.2, except for Fig. The value 
of IZturb/reg is not negligible compared to previous studies -see for example Ref. [29] 
where lZ tur b/reg is low. But the value lZ tur b/reg — 0.2 is still significantly below the 
maximum limit it can reach according to Ref. [27]. 




In the following, we discuss the dependence of the probability to detect the source 
on physical parameters of the source and the magnetic field. For the sources which are 
detected, we find their position on the sky and the local deflection power of the magnetic 
field. 

4-1. Detection of the source: dependence on parameters 

In the following, we define "clusters" as at least three events with energies above 10 19 ' 8 eV 
and located in the same region of the sky, according to the condition given by Eqn. (JSJ) . 
For 5000 events above 10 19 eV in the whole visible sky and no anisotropy, the probability 
to have a cluster of background events is rare (~ 1% only), whereas the probability to 
have doublets of background is of the same order than for source doublets. Therefore 
doublets of events and clusters above 10 19,8 eV have to be studied separately. For 
clusters, the method of Section [3] is not useful because it is optimized to reduce the 
background which is already low in this case with the price of loosing some real sources. 
On the contrary, for doublets, the method is necessary. It reduces the probability 
to detect the source too, but it reduces more strongly the probability to detect the 
background. Therefore we separate doublets and clusters on all plots and mostly focus 
on doublets in this section. However, for higher numbers of events in the sky, clusters 
of background events are more frequent and the method of Section [3] can be used for 
clusters in this case too. 

In Figs. EHHl we plot respectively with the magenta dashed-dotted and the black 
dashed lines, the probabilities to have on the visible sky a doublet of events from 
the source or from the background. This first value gives an estimate of how many 
sources of equivalent luminosities would be required so that one of them would emit 
a doublet. The red thick and the blue thin solid lines respectively correspond to the 
probabilities to detect the doublets from the source and from the background with the 
method proposed in Section [3] The ratios between these two sets of lines respectively 
show the efficiencies of source detection and background rejection, when one starts with 
a high energy doublet. The green dotted lines correspond to the probability to have 
either a doublet or a cluster from the source. The total probability to detect the source 
is shown on two plots with very thick green lines. It corresponds to the sum of the 
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Figure 6. Dependence on source parameters. Left (a): source luminosity (proton 
source). Right (b): CR composition (2.6% for source luminosity). Very thick green 
solid line for the total probability to detect the source. E t ^ = 10 196 eV, 5000 events 
above 10 19 eV in the whole visible sky, B reg = 2 [iG and B rms0 = 4 /j,G. 

probability to detect doublets and of the probability to have a cluster, which has to be 
seen as a detection in ~ 99% of cases. 

The results depend on parameters of the source: its luminosity and the composition 
of primaries. Fig. [6|i shows an example of quantitative results for the dependence on 
the source luminosity, in the range 0.9% to 3.5%, for 5000 events above 10 19 eV in the 
whole sky and E t h = 10 19 ' 6 eV. 0.9% corresponds to a mean value of 0.5 events emitted 
by the source above 10 19,8 eV, and 3.5% to 2 -luminosities of ~ (0.5 — 2) • 10 41 erg-s^ 1 
for sources at 100 Mpc. For a luminosity of 2.6%, the probabilities to have a doublet or 
a cluster from the source are similar because the source is bright enough to emit with 
comparable probabilities 2 or "3 or more" events above 10 19 ' 8 eV. Indeed, for such a 
number of events in the sky, it corresponds to a source which emits a mean value of 
about 1.5 events above 10 19 ' 8 eV. 

If the density of UHECR sources is low, we always can consider that their locations 
on the sky are well separated. For a high density of sources, it is not so, especially in 
the centers of galaxy clusters. However, since the angular resolution of detectors is 1-2 
degrees, all sources located in the center of a galaxy cluster can be considered as an 
effective single source. Such centers of galaxy clusters would serve as effective bright 
sources for our study. A very few nearby galaxy clusters are exceptions to this, because 
they are not point-like -see example of the Virgo cluster in Ref. [T3] . 

In Fig. [6b, we investigate the efficiency of the method for sources of light nuclei 
with Z < 6. For Z > 3, the background detection is larger than the source detection. 
4 He photo-disintegrates at 10 Mpc distances. Then, in these conditions, light nuclei can 
only be detected in case of a bright enough source, thanks to clusters of events at high 
energies -see very thick green solid line. In case of experimental statistics much higher 
than 5000 events above 10 19 eV, one could take higher energy thresholds, which would 
improve the detection of light nuclei sources -see below. For low TZturb/reg, this method 
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Figure 7. Dependence on statistics. Left (a): energy threshold of sky maps. Right 
(b): total number of cosmic rays in the sky, above 10 19 eV. 2.6% for source luminosity, 
B reg = 2 fxG and B rms o = 4 /j,G. 5000 events above 10 19 eV in the whole visible sky 
for (a), and E th = 10 196 eV for (b). 

could also be optimized for and applied to bright sources of heavy nuclei in the future. 

With such statistics and values, there is no visible difference between the two source 
spectra (a = 2.2 and a = 2.7). The lower energy tail around the source does not contain 
enough events to distinguish the power law index. 

The results also depend on statistical parameters: the energy threshold of sky 
maps, E th , and the total number of CR detected by the experiment above a given 
energy which we will set to 10 19 eV for this discussion. Fig. [7^ shows the impact of 
the energy threshold E t h for a proton source, with 5000 events in the whole sky above 
10 19 eV and a luminosity equal to 2.6%. The higher the energy threshold, the lower the 
probability to detect the doublets emitted by the source. However, when the energy 
threshold is decreased, the probability to be confused by the background grows faster 
than the probability to detect the source. Giving a general optimum energy threshold 
for source detection is not possible. For example = 10 19 ' 6 eV ~ 4.0 • 10 19 eV (or 
Fj t h = 10 195 eV ~ 3.2 • 10 19 eV) can be typically good compromises between sufficient 
source detection and sufficient background rejection. Fig. [7b is done for the same source, 
with E t h = 10 19 ' 6 eV and for a total number of cosmic rays in the sky above 10 19 eV 
ranging from 2500 to 10000 events, which corresponds to future values for Telescope 
Array. For such a luminosity, the line corresponding to the detection of the source 
doublet grows slower than the line corresponding to the detection of background. This 
is just due to the fact that when the experimental statistics increase, the probability 
to have a doublet becomes less important for such a bright source and it starts to be 
overwhelmed by the probability to have a cluster. The total probability to detect the 
source does not grow slower than the background -see very thick green line. 

For one order larger statistics (10000-100000 events), one could increase both 
and the energy required for the two events of the high energy doublet (currently set 
to 10 19 ' 8 eV), because the statistics at very high energies would be multiplied by ~ 10. 
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Figure 8. Dependence on GMF parameters. Left (a): regular field strength close 
to the Sun, B reg , with a constant ratio B reff /B rms o set to 0.5. Right (b): turbulent 
field strength, B rms o, with B reg = 2 fiG. 2.6% for source luminosity, E f /j = 10 196 eV, 
5000 events above 10 19 eV in the whole sky. 
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Figure 9. Limitations of the correlation coefficient method: dependence on the ratio 
Hturb/reg of the percentages of detected source and background "2+1" -see text. Thin 
red and thick magenta solid lines for detected "2+1" from the source with G m in = 0.95 
and C m i n = 0.999, respectively. Dotted blue line and dashed black line for detected 
"2+1" from background with C m i n = 0.95 and C m i n = 0.999, respectively. 2.6% for 
source luminosity, 5000 events above 10 19 eV in whole the sky, Ey, = 10 19 6 eV and 
B reg = 2 fiG. V ~ 5.3° x 40 EeV. 



Note that this would improve the efficiency of the method for heavier nuclei. 

We also study the dependence on GMF parameters. In Fig. Ek, we plot the 
probability to detect the source versus the regular field strength, with a constant ratio 
between the regular and the turbulent fields strengths set to B reg /B rms0 = 0.5. In Fig.|8b, 
we show the dependence of the probability to detect the source on the turbulent field 
strength B rms0 , with a fixed regular field strength of 2 /jlG. The five internal parameters 
of the method have been optimized for B reg = 2 ^tG and B rms o = 4 //G. They are 
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kept unchanged for the other points of these figures. As expected, for high values of 
the regular and turbulent fields, the results are worse. However, they are still quite 
acceptable even for the highest field strengths which are considered in Fig. El 

With the parameters discussed above, in 60 % of cases there is only one event in 
the energy range 10 19 ' 6 eV < E < 10 19 ' 8 eV. We call such cases "2+1". The fact we 
require at least three events above already allows to discriminate efficiently between 
source and background doublets. In addition to this requirement, the use of the sector 
is another key point to detect the source. Note that for the remaining 40% of cases, 
the probability to have "2+more than 1" from the background is negligible compared 
to the probability to have it from the source. The method can also be optimized for 
stronger magnetic fields by changing the requirement on the minimal number of events 
in the sector. 

For the case of "2+1" events, one can see the limitations of the correlation coefficient 
method for large values of the turbulent GMF. Fig [9] gives the fraction of "2+1" cases 
from the source and the background that are detected through the correlation coefficient 
method, for different ratios TZturb/reg, and for two different values of C min (0.95 and 
0.999). It shows that one can still strongly improve this discrimination between "2+1" 
from the source and "2+1" from the background by using the correlation coefficient 
method only for low values of TZturb/reg -below ~ 0.1. It is still worth for ~ 0.25, but 
useless for much higher values like ~ 0.4. In this later case, for such statistics one can 
simply count events in the sector and require at least one third event between 10 19,6 eV 
and 10 19 ' 8 eV in order to discriminate between source and background. 

The previous results are computed for an isotropic sky, even at the highest energies. 
If we introduce a reasonable anisotropy in the arrival directions of CR above 10 19 ' 8 eV, 
by assuming that all of them are located in a given fraction of the sky, ranging from 
25% to 100%, we see that the results still change in a linear way. For example, if all 
events above 10 19 ' 8 eV are located in 25% of the visible sky, both the probabilities to 
have a background doublet and to detect it are multiplied by 4. 

4-2. Reconstruction of the source position and of the regular GMF deflection power 

With a turbulent GMF strength set to B rms0 = 4 /xG, the method can reconstruct 
the source position at less than 1 degree from the real one in 68% of cases -see the 
distribution with a blue thick solid line in Fig. fTOk . This result is satisfying because the 
experimental angular resolution is close to 1 degree. 

Also, for this turbulent field value, the deflection power can be computed with a 
precision of about ± 25%, in 68% of cases -see blue thick solid line in Fig. ITOb . 68% of 
all reconstructed values are between 4.5 and 7.4° x 10 19,6 eV, and 95% are between 3.3 
and 12.4° x 10 19 6 eV. The distribution has a bigger tail for higher values of T>. 

The precisions on both results are mostly affected by the turbulent GMF strength: 
the distributions computed for B rms0 = fiG and B rms0 = 2 /zG are shown on the 
same figures for comparison. The lower the turbulent component, the more precise the 
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Figure 10. Left (a): distribution of distances, in degrees, between reconstructed 
and real sources. Right (b): distribution of reconstructed deflection powers in 
degrees x 10 19 6 eV. Thin red solid line for the turbulent component off (0 /iG), 
magenta dashed line for B rmS Q = 2 ^G, and thick blue solid line for B rms o = 4 /iG. 

results. With B rms0 = /iG, in most of the cases, the distance between the real and 
the reconstructed source is much lower than one degree -see thin red lines in Figs. [TOa 
and b. Then, the precision of reconstruction would be completely dominated by the 
experimental angular resolution. The distribution for the reconstructed deflection power 
T> has a small width and is peaked around the theoretical value. 

The value of the regular component -for B reg = 1 to 3 fiG- does not affect the 
precision on T>. For low B reg , the results for source reconstruction are more precise. 

For the values we consider, the detector resolution on energies does not really affect, 
statistically, the reconstruction of the source position. If we take AE/E = 25% which is 
the order of commonly admitted values, both for the distance between the reconstructed 
and the real source and for the reconstructed deflection power, no significant difference 
with AE/E = is found. 

No really noticeable impact of the luminosity -between 1.5 and 3.5%- or the total 
number of CR in the sky -between 2500 and 10000 above 10 mo eV- could be detected 
on the results. 

5. Conclusions 

We investigated here the possibility to identify single bright ultra-high energy cosmic 
ray sources on the sky if the regular component of the Galactic Magnetic Field gives 
the major contribution to deflections of UHECR. 

To simulate the background, we took events following the Telescope Array exposure. 
Their energies were distributed according to the cosmic ray spectrum measured by 
HiRes. We considered the cases of 2500 to 100000 events with E > 10 19 eV. 

For the single source, we took a power law acceleration spectrum proportional to 
1/E 22 , with a maximum energy E max = 10 21 eV. Then, we propagated it from the 
source, put at distances equal to 50-100 Mpc, to the observer by taking into account 
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all energy losses as in Ref. [8J. In order to simplify the study, we assumed, following 
Ref. that deflections in the extragalactic magnetic field are small. We deflected the 
cosmic rays produced by the source in the Galactic Magnetic Field. We considered a 
regular component in the disk and the halo similar to the Prouza and Smida model [19], 
but we modified some parameters according to a recent discussion during the Ringberg 
workshop on UHECR and magnetic fields [25] (for details, see Section |2~7L]) . We also 
added the turbulent component both in the disk and in the halo. Finally, we modified 
the signal from the source according to the exposure of Telescope Array. 

We combined the signals from the source and the background in a common data 
set and tried to recover the source from these data, using the method presented in 
Section [3j As we showed in Fig. [H in order to discriminate between the source and the 
background when the contribution of the turbulent GMF is not negligible, one needs at 
least 2 events from the source at the highest energies, E > 60 EeV. Then, one can look 
for the low energy tail, as discussed in Section [3j 

We studied in Section H] the dependence of these results on the main parameters of 
the source and of the Galactic Magnetic Field. 

Finally, we reconstructed the position of the source and the deflection power of the 
magnetic field in Section 14.21 The main result is that the position of the source can 
be reconstructed within one degree from the real position in 68% of cases. The regular 
GMF deflection power is known with a 25% precision level for a turbulent field strength 
equal to 4 /iG close to the Sun. For lower values of the turbulent magnetic field, the 
precision is much better. 
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